Fractal spin structures as origin of 1/f magnetic noise in superconducting circuits 
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We analyze recent data on the complex inductance of dc SQUIDs that show 1/f inductance noise 
highly correlated with conventional 1/f flux noise. We argue that these data imply a formation 
of long range order in fractal spin structures. We show that these structures appear naturally in 
a random system of spins with wide distribution of spin-spin interactions. We perform numerical 
simulations on the simplest model of this type and show that it exhibits 1/./^^'' magnetization noise 
with small exponent and reproduces the correlated behavior observed experimentally. 
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Despite recent progress, the origin of low-temperature 
(T < IK) 1/f flux noise in superconducting circuits re- 
mains an open question for nearly 30 years 1] . A progress 
in this problem would be important for a number of ap- 
plications: noise suppression would improve the sensi- 
tivity of SQUID magnetometers [l|, [2| and eliminate the 
dominant source of decoherence for flux Q and phase [1] 
qubits. The problem has two parts: the microscopic ori- 
gin of the degrees of freedom that produce the noise and 
the origin of the interaction responsible for their dynam- 
ics. 

It is now accepted that the noise is due to the electron 

fins localized at the surfaces of the superconductors 0- 
The surface localization is implied by the weak depen- 
dence of the noise amplitude on the area of SQUIDs: vari- 
ation of the loop area over 5 orders of magnitudes does 
not lead to a significant systematic change in the noise 
amplitude [1, 5]. The approximate temperature indepen- 
dence of the noise in the range 20mK ^ T < 500mK sug- 
gests that the magnetic system responsible for the noise 
must be characterized by very low energy scales pointing 
to electron or nuclear spins. The nuclear spins can be 
ruled out because 1/f flux noise persists up to lOMHz 
which is much larger than the typical energy scale as- 
sociated with nuclear spins / < IkHz Q. Recent data 
[ij, 9] confirm that the flux noise is due to electron spins 
situated close to the surface of superconductors and show 
that these spins are characterized by the surface density 
ris ~ 5 X lO^^m"^. Further support for this conclusion 
is provided by density functional calculations [lo| which 
found that localized electronic levels are likely to occur 
at strongly disordered metal-insulator (MI) interfaces. 

The dynamics of the electron spins that generate the 
noise remains poorly understood. However, very recently, 
an experimental work by S. Sendelbach et al. |11| discov- 
ered a highly unusual feature of the magnetic noise in dc- 
SQUIDs at millikelvin temperatures that severely limits 
possible mechanisms. It was found that the noise in the 
SQUID inductance has also 1/f power spectrum and that 
the inductance fluctuations are highly correlated with the 
usual 1/f flux noise. Moreover, the correlation between 
the two noises grows as the temperature is decreased and 



becomes of the order of unity below T ~ lOOmK. This 
suggests a common underlying mechanism producing the 
noise in both inductance and flux. Because inductance 
and flux are even and odd under time inversion oper- 
ation respectively their cross-correlation implies broken 
time inversion symmetry and the appearance of a long 
range magnetic order. 

The formation of magnetic order is difficult to reconcile 
with the temperature independence of the flux noise be- 
cause normally the former implies local fields of the order 
of the transition temperature, Tc, that suppress individ- 
ual spin dynamics at T ^ Tc. The collective modes are 
expected to be suppressed by anisotropy at the temper- 
atures in mK range so they cannot lead to temperature 
independent flux noise. Similarly, the transitions between 
different metastable states in ordered magnets involve do- 
main wall motion which disappears at low temperatures. 
Finally, the spin glass formation has to be ruled out be- 
cause one does not expect correlations between magneti- 
zation and susceptibility in spi n g lasses, this expectation 
was confirmed by simulations |12l |. 

In this paper we show that the experimental results fllj 
are reproduced in a highly disordered spin model where 
interaction between spins is characterized by a very broad 
distribution of the coupling strengths. In this situation 
the long range order is due to formation of a large fractal 
cluster of strongly coupled spins which spans the whole 
system. This ordered cluster includes only a fraction of 
spins leaving many smaller clusters and isolated spins 
free to fiuctuate. We formulate the simplest model of 
this type and provide numerical results that agree with 
the existing data. In more detail, we consider a system 
of Ising spins distributed randomly in a 2D plane. Spin- 
spin interactions are assumed ferromagnetic and decaying 
exponentially with distance which provides a broad dis- 
tribution of couplings between them. We find that mag- 
netization dynamics simulated by single spin flip Monte 
Carlo algorithm produces noise with 1/f power spectra. 
We computed the linear response, x(0; of the system 
to a magnetic field of the low frequency, w, and fiuctua- 
tions of the magnetization of the system, M{t), at time 
scales much longer than the period of the external field. 
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tuj ^ 1. These quantities mimic the properties measured 
in the experiments [ill, E^l • Our main result is a large 
cross-correlation between the noise in the susceptibility 
and magnetization, as illustrated by Fig. [1] in a wide 
range of temperatures in the ferromagnetic phase. 

The model provides a new insight into the nature of 
spin-spin interactions which is the most likely mechanism 
generating the dynamics of the spins. First, we note that 
both phonons and nuclear spins can be ruled out as a 
source of electron spin dynamics due to their weak cou- 
pling to the electron spins ^5, 6]. Alternative sources of 
the dynamics such as hopping of electrons between traps 
with random spin orientations or interaction with tun- 
neling two level systems (TLS) [l3| are difficult to rec- 
oncile with experimental data. The former requires very 
high concentration of single-electron excitations while the 
latter implies high concentration of thermally activated 
TLSs . This leaves only electron spin interaction with 
each other. 

The typical spin-spin distance inferred from the re- 
ported densities is > Inm, which implies an exponen- 
tially small direct spin exchange. However, proximity to 
the metal makes RKKY interaction possible This 
interaction can be especially large for the electrons lo- 
calized on metal-insulator boundary levels that remain 
singly occupied due to Coulomb repulsion [l^, while 
electrons localized deeper in the insulator are expected 
to have exponentially weaker coupling to the conduction 
electrons. The exponential dependence of the coupling 
to conduction electrons implies a broad distribution of 
spin-spin couplings. 

Theoretical analysis of a system of spins interacting via 
such random RKKY mechanism under assumption of suf- 
ficiently weak interactions (such that the system does not 
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Figure 1: Cross correlator of noise power spectra of the mag- 
netization and susceptibility for T = 0.05Jo and a particular 
disorder configuration (circles); T = 0.05 Jo averaged over dis- 
order (squares); T = O.lJo averaged over disorder (triangles). 
Inset shows temperature dependence of the normalized dis- 
order averaged correlator of magnetization and susceptibility 
fluctuations (Mx"). 



form an ordered state) showed that the system can gener- 
ate temperature independent 1/f noise. Relatively high 
frequency part of the noise spectra [1| is generated by the 
diffusion of local magnetization fd\ , whereas the very low 
frequency part [4, TSj is generated by fluctuations of rare 
pairs of spins [16l |. 

Temperature independence of the 1/ f noise and corre- 
lations between flux and inductance [ll| imply unusual 
magnetic long range order. In a conventional magnet 
the noise is due to domain wall jumps. Each such jump 
induces changes in the magnetization and susceptibility 
proportional to the changes in the domain volume SV 
and domain wall area 6S respectively. This means that 
in conventional magnets where S <^ V the fluctuations 
in the response to an external field is much smaller then 
fluctuations in the quasistatic magnetization in contrast 
to the data. 

Large fluctuations of the response and as well as high 
correlations between the response and magnetization are 
possible in highly disordered magnets with wide distribu- 
tions of spin-spin couplings where correlated spins form 
fractal clusters with V ^ S. The simplest model that 
contains this physics is characterized by the Hamiltonian 

H = - ^^JijSiSj, Jy = Joexp{~rij/a), (1) 

where i numbers sites of the random lattice, a is a decay 
length of interactions, Jq determines the annealing tem- 
perature of the system and spins are classical s = ±1. 
Three dimensional version of this model was used earlier 
to analyze the formation of l ong range order in magneti- 
cally doped semiconductors (l7l ITsj. 

At small densities na^ <C 1 the system ([1]) undergoes fi- 
nite temperature ferromagnetic transition which is driven 
either by the temperature T or by the lattice site density 
n. This transition is in the universality class of perco- 
lation transition. It can be understood in terms of the 
"circle packing" problem in 2D in which one considers 
circles of radius, r{T)/2 = a/2lnjQ/T, drawn around 
each spin site. Overlapping circles correspond to spins 
separated by Vij < r{T) which are therefore strongly in- 
teracting and, thus, aligned. In contrast, spins separated 
by rij > r(T) are effectively independent. Dimensionless 
parameter B{T) = i:r'^{T)n controls the thermodynamic 
transition. At high temperatures B{T) <^ 1 strongly cou- 
pled spins are rare and form only small clusters. As the 
temperature decreases, B(T) increases, both the num- 
ber and sizes of the clusters grow and at some threshold 
value B = Be an infinite cluster of strongly coupled spins 
forms. This corresponds to the ferromagnetic ordering at 

the critical temperature, — Jocxp {^—^J Be/ (iTna^)^ . 

The threshold value for the "circle packing" problem is 
Be ~ 4.5 [l^. In the vicinity of the peroclating tran- 
sition the cluster size distribution is given by the power 
law: n{s) ~ exp(— Cs"^/^), B > B^ where C is a con- 
stant and r is a characteristic exponent [l9| . The scaling 
behavior is realized for \B — Bc\ <C Be which translates 



3 




Figure 2: Cluster size distribution n{s) at \Bc — B\/ Be = 0.2. 
Up to sizes s > 100 the data fits a power law n(s) ~ with 
r ~ 2. 
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Figure 3: Normalized power spectra of the magnetization 
noise at T = 0.05Jo of the lattice with A'' = 1024 sites. Lin- 
ear fit yields a = 0.93 for the noise exponent. Inset shows 
Smx"U) at T = 0.07 Jo, a = 0.72. 



into a wide range of temperatures: ln(Tc/r) <C Xj \/no? . 

We analyzed numerically the cluster size distribution in 
the model ([T]) in the parameter range < \Bc — B\l ^ 
0.2. This translates into |lnTc/T| < 0.5 for the model 
characterized by realistic values of In Jq/Tc ~ 5. We 
found a broad power law distribution of cluster sizes up 
to s > 100 at |Bc - B\IBc = 0.2, see Fig. [21 indicating 
that in this whole temperature range the spin system 
separates into many large clusters weakly coupled to each 
other and to the infinite percolation cluster. 

Formation of large clusters of spins that flip as whole 
objects under thermal excitation affect the magnetization 
dynamics of the system for two reasons. First, the con- 
tribution to the noise coming from large clusters scales as 
n(s)s^ so their relatively small number is compensated by 
the additional factor of s. Second, the broad distribution 
of the cluster sizes translates into wide distribution of re- 
laxation rates F associated with cluster dynamics which 
results in 1// power spectra of the magnetization noise. 

In order to check this conjecture we have simulated 
the dynamics of the spin system ([1]) in the critical regime 
using single spin flip Monte Carlo dynamics satisfying 
the detailed balance condition. We focus on the fluctu- 
ations of magnetization in the regime below the freezing 
temperature, Tc, at which the infinite cluster is formed. 
To minimize the transient regime we choose for the initial 
configuration that has spins on the infinite cluster aligned 
in one direction whereas all others are distributed ran- 
domly. Very long runs up to 3 x 10"'^" Monte Carlo steps 
(MCS) were performed to analyze long time dynamics in 
the steady state characterized by finite average magne- 
tization. We observe 1//" power spectra of the noise in 
magnetization of the system, Fig[3l with the noise expo- 
nent 0.8 < a < 1.2. 

In order to simulate the susceptibility noise measure- 
ment pT|. we subjected the spin system to a time- 
dependent magnetic field B = Bq cos 27rwt and sim- 
ulated the lock-in detection by applying digital low- 
pass filter [2^ to the fluctuation of magnetization M (t) 
and its Fourier component M{t) * cos(27rwt -I- (p), where 
(j) = 7r/2 corresponds to the out-of-phase response 
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We take the field period to be relatively long, 
= 3 X 10^ MCS, for the system of TV = 1024 
The resulting cross-correlator of the power 
SMx"if)/iSM{f)S^>'{f)y/^ where SmU) - 

Jo 



dte~^-^f^'M{0)M{t), is shown in Fig. [T] 
Temperature dependence of the cross-correlator of the 
noise amplitudes {Mx") averaged over disorder is shown 
in the inset of Fig.[TJ power spectra of the cross-correlator 
is shown in the inset of Fig. [31 

The noise spectrum and correlations observed in nu- 
merical simulations can be understood using the follow- 
ing qualitative picture. Long time scale dynamics respon- 
sible for the correlated noise in susceptibility and mag- 
netization is generated by flipping parts of infinite perco- 
lation cluster and large clusters at its boundary. These 
flips are due to thermal activation, so the rate of the clus- 
ter to flip is F ~ exp(— y/T), where V is the maximal 
energy barrier encountered on the optimal spin flip path 
that flips the whole cluster. 

Fractal geometry of clusters results in a logarithmic 
dependence of V^(s) on the cluster size. This unusual de- 
pendence was first seen in numerical work |23] ; it can be 
understood theoretically by assuming that clusters that 
dominate the dynamics at relatively high frequencies are 
essentially random graphs 2J|. Locally, any such struc- 



ture can be viewed as a tree growing from the 'origin' that 
is connected by to -I- 1 bonds, of strength J, of which m 
connect to to sites of the first level that are connected by 
TO bonds each to another sites of the second level and 
so on (see Fig. [31 for m = 3). The size of the cluster is 
1 + (m^+^ — 1) /(m — 1) where L is the number of levels. 
Denote hy Va — Vb — Vc = V the energy barrier to flip 
each of the subtrees A, B and C in Fig. [4l An optimal 
path to flip the whole tree involves flipping all subtrees 
in sequence which means that the maximum barrier will 
he V + J. This barrier is encountered when one of the 
subtrees, say A, have been flipped which requires keeping 
a bond from "0" to A in the high energy state. Subse- 
quent flipping of B will therefore require V + J. This is 
the maximum barrier encountered because with the next 
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Figure 4: Sketch of the fractal spin structure illustrating the 
logarithmic dependence of the barrier with the cluster size. 

step site "0" can be flipped reducing the barrier for C 
to V -\- J . In the case of m > 3 the maximum barrier 
will be « V + m/2J . This means that the energy barrier 
grows linearly with the number of levels whereas the size 
of the cluster is exponential. Therefore the resulting de- 
pendence of the energy barrier on the size of the cluster 
is logarithmic. One expects the logarithmic barrier scal- 
ing with size to remain valid also in the case of a random 
tree characterized by m and J fluctuating around average 
values (m) and (J), which is a more accurate model. In 
the case of percolation on the diluted Bethe lattice with 
fixed bond strength J = (J) 

y(,)«^ln,, (2) 

where D is the fractal dimension of the tree and the pa- 
rameter Z is restricted by l/ln2 < Z < 2/ln2 21]. A 
very similar argument based on the 'links' and 'blobs' 
structure of the cluster !22] supports the validity of scal- 
ing ([2]) for the whole cluster [2j| . We expect this logarith- 
mic dependence of the energy barriers on the cluster size 
should hold at least approximately for the system of spins 
interacting via exponential interaction ([1]) at percolation 
threshold with a similar value of Z . 

The result ([2|) translates into the power law spectra 



of the magnetization noise. At equilibrium the average 
autocorrelator of the magnetic moment Mi (t) of the i-th 
cluster is given by Mi{0)Mi{t) ^ e^'"**. Treating clusters 
as independent and approximating the cluster distribu- 
tion by n{V, s) ~ S{V — V{s))s~'^ , we get 

At percolation threshold the values of the parameters t = 
187/91 and D = 91/48 are universal 

Out-of-phase response to the periodic magnetic field 
with relatively high frequency lo is dominated by small 
clusters with F ~ w that do not equilibrate during one 
field period. These clusters experience a local quasistatic 
field created by all neighboring clusters that fluctuate 
with frequencies f ^ uj. When one of these clusters 
flips the local field acting on the small cluster changes 
which results in the change in its response to the exter- 
nal field at frequency w. This leads to the susceptibility 
fiuctuations which are correlated with the magnetization 
fiuctuations as observed in numerical simulations and in 
the experiment. 

In conclusion, we have shown that the main results 
of recent experiments (ll| on the flux noise are repro- 
duced by the simplified model of a disordered ferromag- 
net with a very wide distribution of couplings between 
spins. In the simulations we assumed single spin fiip dy- 
namics which does not conserve local spin. Understand- 
ing the mechanism producing such dynamics or general- 
izing these results to other types of dynamics will be a 
subject of future work. 
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